Polyelectrolytes in the presence of multivalent ions: gelation versus segregation 
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We analyze solutions of strongly charged chains bridged 
by linkers such as multivalent ions. The gelation induced by 
the strong short range electrostatic attractions is dramati- 
cally suppressed by the long range electrostatic correlations 
due to the charge along the uncrosslinked monomers and ions. 
A modified Debye-Hiickel approach of crosslinked clusters of 
charged chains is used to determined the mean field gelation 
transition self-consistently. Highly dilute polyelectrolyte solu- 
tions tend to segregate macroscopically. Semidilute solutions 
can form gels if the Bjerrum length Ib and the distance be- 
tween neighboring charged monomers along the chain b are 
both greater than the ion size a. 



Linear polyelectrolytes are ubiquitous in biology given 
that nucleic acids and most proteins are charged. More- 
over, they have important technological applications as 
gelling and drag reduction agents. Long linear polyelec- 
trolytes are typically water soluble in low ionic strength 
monovalent salt solutions due to a net repulsion between 
the charged monomers. Linking agents such as multiva- 
lent ions of valence z > 1, however, modify the stability 
of polyelectrolytes aqueous solutions [1-5]. The stabil- 
ity of biopolymers in various ionic media has important 
implications in biotcchnological processes [6]. 

The solubility of linear polyelectrolyte in multivalent 
salts has been extensively studied in dilute solutions of 
continuously charged flexible chain [7] and related sys- 
tems [8]. Only few studies, however, have discussed the 
stability of the chains in semidilute solutions [2,3]. In par- 
ticular, the strong binding between small (dehydrated) 
multivalent metallic ions and charged (0~) groups in dif- 
ferent chains opens the possibility of a gelation transition 
in semidilute solutions [9,10]. The competition of gela- 
tion and phase separation has been recently studied in 
neutral chains with linkers [11,12]. The charge along the 
chains, however, will strongly modify the phase diagram. 
Here we describe the gelation by multivalent ions includ- 
ing long and short range correlation self-consistently in 
the analysis. 

Our model treats the strong short range correlations 
driven by the multivalent cations explicitly using bridges 
or crosslinks that change the structure function of the 
system and therefore the long range electrostatic corre- 
lations. The long range correlations are accounted via 
a modified mean field approach known as the Random 
Phase Approximation (RPA) on the linked system. RPA 



describes properly long range correlations in polyelec- 
trolyte solutions without short range correlations [13]. 
Therefore, by using as a reference state of RPA a strongly 
correlated system of chains linked by multivalent ions, we 
can in principle determine self-consistently the gelation 
transition in the system of charged chains. We consider 
only flexible charged chains to avoid difficulties found in 
the gel structures of charged rigid rods [14] due to the 
strong oricntational dependent electrostatic interactions 
[15,16]. 

Consider a polymer solutions where every chain 
monomer carries a negative charge. The concentration 
of monomers is p m , and of monovalent counterions is 
Pm = Pm, and we assume for simplicity that all of them 
are completely dissociated from the chains. We denote 
the concentrations of multivalent ions and their coun- 
terions as pi and p|, respectively (p? — zpi). Let us 
assume that z — 2. The interaction between the mul- 
tivalent ions and the chains, in principle, leads to the 
formation of monocomplexes (mc) and dicomplexes (dc) . 
A mc is formed when an ion binds onto one negative 
site of the polyelectrolyte. Formation of a dc brings 
two negative sites together. The competition between 
the formation of mc and dc is described by the process 
2mc < — ► dc + free ion. The equilibrium constant p of 
this process can be estimated as [3] p = K 2 expels) /K 2 , 
where Ib — e 2 /eksT is the Bjerrum length and the in- 
verse screening length k is given by k 2 — AttIb^I z aPa 
with the sum taken over all dissociated charges. K\ and 
Ki are the statistical weights of mc and dc respectively. 
Assuming that all ions and the monomers of the chains 
have the same size a we can estimate these weights fol- 
lowing Fuoss [17] 
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As shown in Ref. [3] dicomplexation is favorable if p > 
1/12. Since in our work we consider strongly charged 
chains regimes with p ~ 10, we neglect the possibility of 
forming monocomplexes. 

Dicomplexation processes may have two different sce- 
narios. If the two sites that form a dicomplex belong 
to the same chain then this chain collapses onto itself. 
If these sites, however, are from different chains then a 
crosslink is formed. The monomolecular collapse is ob- 
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served when the multivalent ions nearly neutralize the 
chains at [7] 



coll _ Prn 

Pl ~ z 



(3) 



Here we consider the formation of branched structures 
by crosslinks at p% < p\ . 

We write the free energy of the system in the form 



F = F* + F iei + F d 



(4) 



F* accounts for excluded volume interactions and can 
include small effective short range monomer attractions 
(x < 2) to account for slightly hydrophobic chain back- 
bones. It can be approximated as 
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TV 



= (l-^)log(l-0)-x^ 
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Here T is the temperature in units of ks, V is the vol- 
ume, <f) m = a 3 p m , and cj> is the total volume fraction of 
all components = a 3 (p m + p c m + pi + pf). F re f cor- 
responds to free energy of the system without any long 
range electrostatic interactions, the reference system, 



fief 

TV 



^77 = Pc ln 

{C} 



PC'SC 

ewe 



(6) 



The sum is taken over all possible structures C of clus- 
ters formed due to crosslinking ({C} also includes trivial 
clusters such as free ions and counterions), and pc, wq 
and sc represent concentration, statistical weight and 
symmetry index of a cluster C, respectively. 

The last term F e i in Eq. (4) takes into account all long 
range electrostatic effects and can be written in the form 
[18] 
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where the matrix of correlation functions gk = \\g^\\ has 
the following form 



(8) 



Here ff(C)J? is the scattering factor of a cluster charac- 
terized by a structure C 



»(tf)i?' = E( e<k( 



(9) 



with to and n running over all i-type and j-type units in 
C cluster respectively. 

We calculate free energy of the reference system as- 
suming that all complex clusters have only tree-like ar- 
chitectures [3,19,20] 



F rof = fi> c d f + F r c e r b + fie f oss (io) 

F r ^ f is associated with all translational entropies, 

TV = ~W ln Pm + P " 1 ln P " 1 + Pl ln Pl + P ' ln ^ ^ 

F r c ™ b is the term obtained from the number of possibili- 
ties to choose the monomers and the ions that participate 
in formation of crosslinks, 
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= Pm (T m inr m + (i - r m ) in(i - r m )) + 

PiiTilnTi + (1-1^)111(1 -TO) (12) 



where T m and Ti are the fractions of the monomers and 
the ions respectively that belong to the crosslinks. Be- 
cause we neglect the formation of monocomlexes these 
fractions are related by p m T m = zpiTi. The last term in 
Eq. (10) corresponds to the crosslinking free energy, 
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We use Eqs. (10)-(13) for the reference free energy to 
determine the instability of the system due to the forma- 
tion of the infinitely large network or gel. The gelation 
line is identified by the divergence of the cluster's weight 
average N w . We note that N w — ► oo corresponds to a 
singularity of the second derivatives of the reference free 
energy [19] 
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Here the set {p a } includes p u p m and p^ oss = p m r m . 
Representing fi e f only as a function of these three con- 
centrations and evaluating the determinant (14) we get 
the critical value that corresponds to the formation 
of the network 
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(z-l)(N-l) 



(15) 



The last term of the free energy F e i given by Eq. (7) 
is calculated as follows. First we evaluate the 5x5 ma- 
trix of correlation functions gk which has the following 
components: gu — p c m is the concentration of counteri- 
ons dissociated from the polymer chains, 522 = jO»(l — 1\) 
is the concentration of free multivalent ions, 533 = p\ 
is the concentration of counterions dissociated from the 
multivalent ions, and the correlation functions between 
the connected units are denoted by g 44 , g i5 — g 5i and 
g 55 , where 4 and 5 stand for non-crosslinked monomers 
and crosslinked aggregates respectively. In the present 
work we neglect dipole-dipole and dipole-charge interac- 
tions, therefore, 345 and g 55 have no effect on the free en- 
ergy of the system. We calculate the monomer-monomer 
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correlations 544 using the diagrammatic technique de- 
scribed in Ref. [21]. By introducing the concentration 
t of all branching structures that could be attached to 
a chain monomer, and the concentration T of all non- 
associated monomers, we obtain the total concentration 
of monomers p m = tT, and the fraction of associated 
monomers r m = 1 — 1/t. By finding the correlation 
function between two monomers that belong to the same 
chain = (T/t)(g(k) — 1) and taking into account 

all possible arrangements of linear chains into tree-like 
structures we get 



544O) = jO m (i - r m ) 



l + (l-zr m )(g(fc)-l) 
l-T m (z-l)(g(k)-l) 



(16) 



where g(k) is the structure factor of a chain defined as 

■9( fc ) = 4E( e ~ lk(r, ~ rj) > 



Strongly charged chains in scmidilutc solutions are 
stretched at length scales of few monomers n due to the 
locally unscreened electrostatic repulsions. We assume 
they behave like rods on short length scales, although 
simulations reveal slightly less stretched local conforma- 
tions [22]. On larger scales the chains obey Gaussian 
statistics with N/n segments per chain. For g(k) we use 
the approximated expression 



g(k) - 1 = 
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+ 



N-n 



l + nkb l + Nnk 2 b 2 /12 



(18) 



where b is the average distance between the neighbor 
monomers along the chain, which gives the right limits 
and it is reasonably accurate. Though in ref. [22] it is 
argue that n is nearly <p m independent, we estimated it 
here as [23] 
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to obtain a bound on the gelation in the largest possible 
perturbed chain conformation. We also calculate here 
the gelation for the least perturbed chain conformations, 
Gaussian chains, by putting n = 1 in Eq. (18); ideal 
chain statistics are observed in charged chains in certain 
hydrophobic backbone regimes [24]. 

When the fraction of associated monomers T m reaches 
the critical value given by Eq. (15) then the correla- 
tion function c/44 (k) diverges at k = 0. This divergence 
represents correlations on the infinite scale and, as dis- 
cussed above, corresponds to the network formation. 

To obtain the long range electrostatic contribution to 
the free energy given by Eq. (7) we introduce the follow- 
ing interaction potential between the charges [25] 



U ij (r) =z iZj l B 



1 - e- r / a 
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which allows us to account for short range repulsions due 
to the hard-sphere nature of the charges. 

Taking the Fourier transform of U % ^{r) we write the 
electrostatic free energy (7) in the form 
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where k 2 = p m + zpi + z 2 pi(l — I\) and 



U(k) 



k 2 (l + a 2 k 2 ) 



(21) 
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We obtain the fraction of crosslinked monomers in the 
system T™ ln , for given p m , pi and ls/b by minimizing 
the total free energy with respect to T m . Because the 
correlation function (744(0) diverges at T m — r*„ we find a 
minimum of the free energy in the interval T m G [ 0, r^J. 
If the minimum is reached inside this interval then the 
system contains only finite size clusters. If the minimum 
is reached at T™ 111 = then an infinitely large network 
is formed in the system. By numerically solving 



TZ n (php m ,K(0) = r* m 



(23) 



we obtain the gelation line p*(p m ). For a given p m , if 
Pi < p* the system contains only finite clusters, if pi > p* 
an infinitely large network or gel is formed. The results 
are shown in Figure 1. Above the straight line which 
corresponds to Eq. (3) the chains are collapsed. The 
dash line labelled as 1 is the gelation line when the long 
range electrostatic interactions are not included (i.e., set- 
ting F e i = in Eq. (21)), and it is independent of the 
chain structure function in the mean field model . Line 

2 is the gelation transition obtained for locally stretched 
chains (the structure factor given by Eq. (18)), and line 

3 is the gelation obtained for Gaussian chain statistics at 
all length scales (n = 1). The gelation transition changes 
by orders of magnitude and depends on the chain con- 
formation when long range electrostatic interactions are 
included. 

In Figure 2 we show the gelation lines obtained for dif- 
ferent values of ls/b and b/a. Notice that the ion size a in 
K2 in Eq. (2) is a parameter different than the average 
distance between the neighboring monomers along the 
chain b. It can be varied to account for the unknown size 
of the ions in the crosslinks and the unknown dielectric 
constant of the local medium around the links. For a = b 
the gelation occurs at rather high pi concentrations. This 
may explain why in solutions of strongly charged poly- 
electrolytes, such as DNA and Polystyrene-Sulphonate, 
where both metallic and organic multivalent ions are hy- 
drated around the chains (a > b), z > 3 are required to 
observe precipitation, and the gelation is not observed at 
p* ~ p m j (Nz(z — 1)) even in scmidilutc solutions but at 
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much larger values of pi. Instead in polyclectrolytcs with 
monomers with ionizable OH groups such as acrylate 
groups, since metallic ions are dehydrated when form- 
ing dicomplexes, divalent metallic ions do precipitate the 
chains, and gels are formed in semidilute solutions [10]. 

We investigated the stability of the system to macro- 
scopic phase separation, and found no instabilities in- 
duced by electrostatics below the gelation lines for the 
range of parameters used here ( ls/b ~ 2, \ < 2). Thus 
we conclude that there is no segregation transition that 
competes with the gelation in the regimes studied here. 

In conclusion, we have determined the effect of long 
range electrostatic interactions on gelation of polyelec- 
trolyte chains induced by multivalent ions. The approach 
can be used to describe gelation of charged chains bridged 
by other type of linkers through modification of K2 ■ We 
could also account for variable weights of different types 
of crosslink functionalities m by allowing a K m depen- 
dence. It is of future interest to estimate the effect of 
cyclization of the chains in the gelation theory of poly- 
electrolytes. 
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FIG. 1. Phase diagram for the solution of polyelectrolyte 

chains in the presense of divalent ions. Concentrations are 

given in o 3 units, Isjb = 2, b = a, N = 100. 
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FIG. 2. Gelation lines for different values of parameters 
ls/b and b/a. From top to bottom (ls/b, b/a) — (2,1), 
(3,1), (2,1.5), (3,1.5). Concentrations are given in a 3 units, 
N = 100. 



